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ABSTRACT 

To study the evolution of protoplanetary dust aggregates, we performed experiments with 
up to 2600 collisions between single, highly-porous dust aggregates and a soHd plate. The dust 
aggregates consisted of spherical Si02 grains with 1.5 /im diameter and had an initial volume 
filling factor (the volume fraction of material) of (/jq = 0.15. The aggregates were put onto a 
vibrating baseplate and, thus, performed multiple collisions with the plate at a mean velocity 
of 0.2 m s~^. The dust aggregates were observed by a high-speed camera to measure their size 
which apparently decreased over time as a measure for their compaction. After 1000 collisions the 
volume filling factor was increased by a factor of two, while after ^ 2000 collisions it converged to 
an equilibrium of « 0.36. In few experiments the aggregate fragmented, although the collision 
velocity was well below the canonical fragmentation threshold of ^ 1 m s^^. The compaction of 
the aggregate has an influence on the surface-to-mass ratio and thereby the dynamic behavior 
and relative velocities of dust aggregates in the protoplanetary nebula. Moreover, macroscopic 
material parameters, namely the tensile strength, shear strength, and compressive strength, are 
altered by the compaction of the aggregates, which has an influence on their further collisional 
behavior. The occurrence of fragmentation requires a reassessment of the fragmentation threshold 
velocity. 

Subject headings: Solar Nebula, Planetesimals, Experimental Techniques, Collision Physics, Solar System 
Origin 



1. Introduction 

The formation of planets in the accretion disks 
around young stars starts with the growth of 
(sub-) micrometer-sized dust grains. Embedded in 
the thin gas of the disk, the dust grains collide 
due to a Brownian relative motion and inevitably 



laboratory experiments have shown that for these 
increasing collision velocities, aggregates are re- 
structured and grow to non-fractal but still very 



porous bodies ( 


Dominik and Tielens 


1997 


Blum 


and Wurm 2000 


1 so that growth can continue to 



stick at the small collision velocities (Blum et al. 



2000 Krause and Blum 2004 1. By this mecha 



nism, dust grains can grow to fractal aggregates 
of ~ 100 /im before systematic drift significantly 



increases the collision velocities ( Weidenschilling 
and Cuzzi|1993 ) . Numerical simulations as well as 



larger sizes. 

Once the aggregation has reached millimeter 
sizes, further growth due to sticking collisions be- 
tween similar-sized dust aggregates slows down. 
Different experiments have shown that collisions 
of mm-sized aggregates result in bouncing or frag- 
mentation: Blum and Miinch ( 1993 1 performed 
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collision experiments with ZrSi04 dust aggregates 
with a volume filling factor (the fraction of vol- 
ume filled with material) oi (p — 0.26. For ve- 
locities exceeding ~ 1 m s^^ they found frag- 
mentation as the dominant process, whereas for 
smaller velocities the aggregates bounced. [HeiBel- 



mann et al. ( 2007 1 performed similar experiments 
with highly porous aggregates {4> = 0.15) of 1 - 
5 mm diameter, which collided with a dusty tar- 
get at a velocity of 0.2 m s~^ or with each other at 
0.4 m s~^. In both types of collisions the results 



were dominated by bouncing. Langkowski et al 



(2008) [hereafter paper II], performed impact ex- 
periments of the same highly porous dust aggre- 
gates of 0.2 to 3 mm diameter and different materi- 
als onto equally porous targets of 25 mm diameter. 
For intermediate velocities (0.5 to 2.5 m s~^) and 
projectile sizes (0.5 to 2 mm) they found bounc- 
ing of the aggregate which was even more likely 
if the target was "molded" to a non-flat surface 
(see Fig. Ic in paper II). Smaller projectiles and 
lower velocities led to sticking of the dust aggre- 
gate on the target surface, whereas larger projec- 
tiles and higher velocities resulted in a deep pen- 
etration of the projectiles with no possibility for 
escape (paper II). All those experiments were per- 
formed without the influence of gravity (free fall, 
drop tower or parabolic flight) and in all experi- 
ments the coefficient of restitution e (the velocity 
after the impact divided by the velocity before im- 
pact) was rather small (e < 0.4). 



Many experiments of paper II and HeiBelmann 



et al. ( 2007 ) were performed with the same sample 



material. While in paper II different compositions 
were used (Si02 (monodisperse spheres), irregular 
Si02, and irregular diamond), HeiBelmann et al. 



(2007) focussed on porous aggregates of 1.5 /xm di- 
ameter Si02 monospheres. They were formed by 
the random ballistic deposition mechanism intro- 



duced by Blum and Schrapler (2004) and Blum 



et al. (2006) [hereafter paper I], had a diameter 
of 25 mm, and a volume filling factor of (f> = 0.15. 
The dust samples could be cut or broken into mm- 
sized aggregates to perform the collision experi- 
ments. Although pure silica is not the most abun- 
dant material in protoplanetary nebulae, paper I 
and II showed that the material properties of the 
dust (composition, size distribution, and shape) 
do not significantly alter the experimental out- 
come. Therefore, we regard the material repre- 



sentative according to their mechanical properties 
of aggregates made of a broad collection of refrac- 
tory grains and use the same material in the ex- 
periments presented below. 

The motivation for this work is the explanation 
of the coefficient of restitution in the collisions of 
the described mm-sized Si02 aggregates, which is 
as low as e = 0.2 in aggregate- aggregate as well as 



in aggregate-target collisions (HeiBelmann et al. 



2007 ) . This means that only a few percent of the 
translational energy are conserved, whereas the 
bulk of energy is dissipated in an unknown man- 
ner. Although the aggregates of [HeiBelmann et al.| 
( [200 7') do not show apparent deformation after 
the collisions, the obvious assumption is that the 
energy is consumed by compression. In a micro- 
scopic view, compression results from the rolling, 
breaking and restructuring of inter-grain contacts 



which dissipates energy (e.g. Paszun and Dominik 
2008). In a macroscopic view, compression can 



be described by a compressive strength curve p(^) 



( [Blum and Schrapler|20d4| [Giittler et al.|2009[ pa- 
per I). The dissipated energy in this context is 
AE = p AV, where AV^ is the decrease of a vol- 
ume inside the dust aggregate with internal pres- 
sure p. 

In this work, we describe an experiment, in 
which a dust aggregate with ^ 2 mm diameter is 
placed on a vibrating baseplate and is thus forced 
to perform multiple collisions with this plate. Al- 
though this is not a zero gravity experiment, grav- 
ity is not important for the individual collisions as 
HeiBelmann et al. (2007) have shown that even 



under microgravity conditions sticking never oc- 
curs at the relevant velocities of ~ 0.2 m s~^ and 
that the collision time is too short for substantial 
gravitational influence (see the image sequence in 
Figs. 4 and 5 of [HeiBelmann et al.| ( p007| ) . The 
velocities are the same as those between aggre- 
gate and target in the experiments of |HeiBelmann[ 
et al. (2007). Performing multiple collisions, the 



cumulative compaction is larger than in a single 
collision and can, thus, be assessed. 

In section [2] we describe the experimental setup 
for the measurement of the aggregate compaction 
in multiple collisions. The results are presented in 
section |3] and the relevance and consequences for 
dust aggregates in the protoplanetary nebula are 
discussed in section IH 
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2. Experimental Setup 

As a starting material for a mm-sized, highly 
porous protoplanetary dust aggregate we chose the 
well-defined dust aggregates introduced and char- 
acterized by Blum and Schrapler (2004). They 



consist of 1.5 /im diameter monodisperse Si02 
spheres, possess a volume filling factor of (j)o = 
0.15 and were formed by random ballistic depo- 
sition um and Schraplerl 120041) . From these 
25 mm aggregates we cut out cubes of ^--^2.5 mm 
diameter with a compaction of the aggregate rim 
to a maximum of (jj — 



0.16 (jHeiBel mann et aL] 
Each small aggregate was weighed (mass 



2007) 



jumping particle 

%.•-...#. 



— V- 



vibrating baseplate 

Fig. 1. — Experimental setup: The dust aggre- 
gate in the box is continuously colliding with the 
vibrating baseplate , while it is observed by a high- 
speed camera. The aluminum baseplate vibrates 
with a frequency of 100 Hz with a peak-to-peak 
amplitude of 0.4 mm. 

mo) and put into a box with plexiglass walls and 
an aluminum baseplate of 40 x 40 mm^ size (Fig. 
[ij . The aggregate was observed with a high-speed 
camera at a frame rate of 380 frames per second 
in back-light illumination while the box was vi- 
brated in the vertical direction at a frequency of 
100 Hz with a peak-to-peak amplitude of 0.4 mm 
for different durations (10 to 80 s). This led to a 
jumping motion of the aggregate and 200 to 2600 
bouncing collisions, after which the aggregate (fi- 
nal mass mend) Eind the eroded material in the box 
(total mass of the debris mer) were again weighed. 

In each image, the cross-sectional area and the 
position of the aggregate were measured. The po- 
sition measurements yield the number of collisions, 
which are underestimated by ~20 %, because col- 
lisions near the edge of the box were not illumi- 




Fig. 2. — Sequence of a bouncing aggregate with 
~ 150 collisions between two successive images. 
The aggregate is rounded first and is then clearly 
getting smaller. The width of the single frames is 
4.6 mm. The example is taken from experiment 1 
in Table HI 

nated well. From the maximum height between 
two collisions the aggregate's velocity at the time 
of the collision can be determined assuming a free 
fall motion of the aggregate. Although the exper- 
iments were performed in air, the frictional effect 
of the ambient gas is negligible, as the maximum 
friction force Fair — GTTrjRv — 9.5 • 10~* N is much 
smaller than the projectile's weight Fgrav — fng — 
3.9 • 10~^ N. Here m, R, v, and rj are a typical 
aggregate mass (m « 4 mg), the corresponding 
aggregate size {R ~ 1.5 mm), an average aggre- 
gate velocity {v « 0.2 m s~^), and the viscosity of 
air (?7 = 17.2 ^Pa s); g is the gravitational accel- 
eration. 

As the velocity of the baseplate is unknown 
for the exact time of collision, a statistical colli- 
sion velocity distribution is presented in section |3] 
The imaged cross-sectional area of the aggregate 
was converted into a volume by assuming a sphere 
which is a coarse approximation at the onset of 
the experiments but a good approximation after 
^^150 collisions (Fig. [2]). A total of 18 individual 
bouncing experiments was performed (Table [T]) . 

3. Results 

In this section, we present the calculation of the 
distribution of collision velocities. We will then 
quantify the compaction of the dust aggregates 
and give an analytical approximation for practi- 
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Table 1: Experimental results for all performed experiments. In the experiments marked with ^ the aggregate 
fragmented shortly after the given number of collisions. 
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cal use. Moreover, we will present further conse- 
quences of the structural change of the aggregates, 
namely the fragmentation of a dust aggregate at 
small velocities and the development of the coeffi- 
cient of restitution. 

3.1. Collision Velocities 

The maximum height h of the dust aggregate 
between two collisions can be measured from the 
image with the highest aggregate position and it 
determines the aggregate's velocity at the time of 
impact, V = y/2gh. As the impact time is only 
known with an error of 2.6 ms , as the maximum 
plate velocity of 0.13 m s~^ is in the same range as 
the collision velocity, and as the baseplate velocity 
rapidly changes, we make a statistical approach 
to calculate the distribution of collision velocities. 
The probability of a given plate velocity is given 
as 

P{v) dv ^ v ■ {t{v) - t{v + dv)) , (1) 

where v — 100 Hz is the oscillation frequency of 
the plate and dv determines a velocity interval 



around v. t(v) is the inverse velocity cosine func- 
tion 

t{v) = — arccos ( — ^ ) , (2) 
uj \AoUjJ 

where Aq — 0.2 mm is the amplitude of the plate 
and ll! — 2-Kv is the angular frequency. For an 
aggregate with velocity Uag. plate velocities v < 
—Wag. do not lead to a collision while the maxi- 
mum plate velocity Vpi, is the likeliest. Thus, the 
velocity distribution of the plate (Eq. [ij is con- 
volved with a linear collision probability, cropped 
for V < —fag. and shifted by fag.- This is the veloc- 
ity distribution of a single collision. The same pro- 
cedure is performed for each collision and all nor- 
malized distributions summed up yield the overall 
velocity distribution for one experiment . Figure [3] 
shows the mean of the velocity distributions of all 
experiments (solid line) with the standard devia- 
tion (grey shaded area). All velocities are smaller 
than 0.35 m s^^, the median velocity is 0.19 m s~^, 
while 50 % of all collisions are within 0.14 and 
0.23 m s^i. The velocity range for the individual 
experiments given in Table fl] is the 50 % range 
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Fig. 3. — The mean of the impact velocity proba- 
bilities of all experiments (solid line) and the stan- 
dard deviation from this distribution (grey shaded 
area) . 50 % of the collision velocities are between 
the dashed lines, 90 % are between the dotted 
lines. 

from the individual distributions. 

3.2. Compaction of the Dust Aggregate 

The prime objective of the experiments is to 
measure the increase of the averaged volume filling 
factor of the dust aggregate after n collisions, i.e. 



(n) = 



m{n) Vo 



niQ V{n) ' 



(3) 



Here, toq is the initial mass of the aggregate, 
and Vo is the initial volume calculated from mg 
by assuming an initial volume filling factor of 
00 = 0.15 (Blum and Schrapler||2004|). The vol- 



ume V{n) was calctilated from the projectile's 
cross section (7a(n) by assuming a sphere, thus, 
V{n) = |7r~^/^(Ja('^)^^^- The mass m(n) slightly 
decreases due to erosion at the target. However, 
experiments with different durations did not show 
any change in the relative erosion (see Table [ij 
and most of the eroded fragments were sticking 
to the baseplate at the place of the first few colli- 
sions. Thus, assuming that the erosion took place 
in the very first collisions and was caused by the 
preparation of the sample, we take m(n) = mond- 
The decrease of the volume of the dust aggre- 
gate over the course of collisions with the base- 
plate is obvious as presented in the image sequence 
of Fig. [2] From each experiment we calculated 
the volume filling factor (t>{n) from Eq. [s] and 




500 1000 1500 

number of collisions 



2000 



Fig. 4. — Mean increase of the volume filling fac- 
tor from 18 individual experiments (solid line) and 
the standard deviation (grey shaded area). The 
dashed line represents an analytic approximation 
(Eq. |4]) which converges to a filling factor of 
4> = 0.37. The deviation for small collision num- 
bers is determined by systematic errors in the ex- 
periments while the analytic curve is constrained 
to (j){n = 0) = 00- 

took a boxcar average over 100 collisions to reduce 
stochastical scattering from the rotation of the ag- 
gregate. The mean filling factor and the standard 
deviation of all experiments are presented in Fig. 
|4] (solid line and grey shaded area) . 

The filling factor for n < 150 is greater than 
the initial filling factor of (I)q — 0.15, which is 
due to a combination of various systematic errors: 
The scaling is inaccurate if the aggregate is in the 
front or back of the box which results in an er- 
ror of ±3 % in 4>, and the choice of the threshold 



+5 % 
-0 %■ 



In 



for estimating the size yields an error of 
three experiments, the aggregate appeared from 
the smallest side in the beginning which overesti- 
mates 0, and the assumption of a sphere instead 
of the cuboid underestimates the filling factor un- 
til the aggregate is significantly rounded. Due to 
those uncertainties it is reasonable to regard the 
data only for n > 150 and assume (f){n = 0) = 0o- 
We give an analytic approximation which rep- 
resents the filling factor for n > 150 (dashed line. 
Fig. H: 



with 



= 0.365, A4> 



(4) 



and V 



700. Accounting for the systematic error in the 
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collision number n, which is underestimated by 
20 % (section [2]), we take v = 850 for practical 
use. For later application, we can calculate the 
volume of the aggregate (from Eqs. [s] and |4]) as a 
function of the collision number 



Vin) 



00 ■ Vq 

- A(j) ■ e^"/" ' 



(5) 



where we ignore the mass loss in the first collisions 
and take m(n) = mg. 

In a very simple model, we assume that the 
compression is the cause for the loss of kinetic en- 
ergy, thus 



If we assume that this relation holds for every col- 
lision, we can calculate the volume fraction of the 
compressed material for the equilibrium situation 
71 — !■ 00 to be: 



/c = 



/ AFc dn 
V{n —f 00) 

■ f^Vdr 



(11) 
(12) 

(13) 
(14) 



1-e^ 



m ■ V — p ■ AV, 



(6) 



where AV for the first collision can be calculated 
by deriving Eq. [5] for n — and An = 1: 



AV = 



(7) 



Thus, we can calculate the pressure in the aggre- 
gate (first collision) as 

^ ^ (l^s^WnVM^ ^ (8) 

2 ■ Vb • (0max - M 



For an aggregate with 



4.24 mg, Vo 



14.1 mm'^, V = 0.2 m s , and e — 0.2 this yields 
a pressure of 3424 Pa. Using the compressive 



strength curve proposed by Giittler et al. (20091 



we can calculate the volume filling factor 
compressed volume AVc to be 



exp 



5P-Ig p„ 
A 



1 



in the 



(9) 



For omindirectional dynamic compression, |Giittler| 
et al. (20091 developed the empirical parame- 
ters 01 = 0.12, 02 = 0.58, A = 0.58, and 
Pni — 1300 Pa, which yields a volume filling factor 
(f>c = 0.43, which is slightly higher than the end 
compression ^max = 0.365 of the aggregate. 

The relation between the volume decrease of 
the agglomerate AV and the volume AVc with 
compressed material is determined by the mass 
balance in the volume {AV + AVc): 



AVc 
AV 



(10) 



Inserting the given values, we learn that the vol- 
ume of the compressed material for n ^ cx) is 77 % 
of the end volume of the agglomerate. Thus, the 
aggregate is inhomogeneously compacted and has 
an unaltered core of 61 % in radius. 

3.3. Influence of the Compaction on the 
Mechanical Properties 

The change of the volume filling factor must 
clearly have an influence on the mechanical prop- 
erties of the dust aggregate. One related finding 
is that in four experiments the aggregate frag- 
mented during the succession of impacts (Fig. |5|. 
Those experiments are marked in Table [T] and do 
not show a clear systematic difference to the other 
experiments in which the aggregate did not frag- 
ment . The number of collisions and the filling fac- 
tors in those experiments are rather high - except 
for experiment 16. A possible explanation is a de- 




Fig. 5. — The foiu' cases in which the aggregate 
fragmented (from The fragmentation is rather gen- 
tle, such that the aggregates only break into few 
The width of a single frame is 7.7 mm. 

crease of the critical fragmentation velocity with 
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increasing volume filling factor or, at least, an in- 
creasing breakup probability. Figure [5] shows very 
gentle fragmentation in contrast to the broad size 
distribution of fragments found at higher veloci- 
ties fBlum an d Miinch|1993 1. This breakup is usu- 
ally found when fragmentation occurs at velocities 
near the fragmentation threshold. However, as the 
dependence on the number of collisions does not 
show a clear tendency, a second explanation is a 
general breakup probability for which we can give 
a rough estimate for mm-sized dust aggregates in 
low- velocity collisions as: 

4 ± 2 

Pfr^g = — = (1.8 ± 0.9) • 10-* per collision, 

(15) 

where ^ ncoii — 22650 is the total number of col- 
lisions in all experiments (Table [l]). 




500 1000 1500 

number of collisions 



2000 



Fig. 6. — The jump height of the aggregate (mean 
of all experiments) over the number of collisions. 
The averaging was taken for as many experiments 
as possible: While for the first collisions all exper- 
iments could be utilized, for collision numbers of 
more than 2000 only few experiments were avail- 
able. For n > 150 the jump height only slightly 
increases with the number of collisions (also mean- 
ing compactness), which unexpectedly indicates a 
constant coefficient of restitution. 

Another influence of the structural change 
might be expected for the coefficient of restitu- 
tion: The coefficient of restitution is a measure 
for the dissipation of energy and from section 3.2 
we know that compaction is a plausible dissipa- 
tion mechanism. However, as the compaction is 
not constant over time, the coefficient of restitu- 
tion is not expected to be either. As the mean 



baseplate velocity does not change with time, a 
change in the coefficient of restitution would re- 
sult in a variation of the maximum height of the 
aggregate. Figure [6] shows that the jump height 
of the aggregate does only slightly increase with 
time for n > 150. A linear fit of all heights for 
200 < n < 2000 yields a mean slope of 40 /im / 
1000 collisions. The increasing height in the very 
beginning might be due to structural changes but 
concurrently the aggregate always collides with 
a broad side in these first collisions, which must 
have a substantial but unknown infiuence on the 
height. 

4. Discussion 

In this section we discuss the relevance for the 
protoplanetary nebula, namely, scaling the results 
for different sizes and velocities and estimating 
whether multiple bouncing collisions can occur in 
a reasonable timescale. Furthermore, we discuss 
the consequences of the aggregate compaction for 
their further evolution. 

4.1. Collision Model 

We will develop a scaling relation for the ag- 
gregate compaction in size and collision velocity. 
From the momentum balance of the colliding ag- 
gregate, we can give the pressure in the aggregate 
as 

(1 + e) ■ m ■ V 



P = ' (16) 

T • A 

T-4-7r-i?2 ' 

where (1 -I- e) ■ m ■ v is the change of momentum 
of the colliding aggregate, taking place within the 
collision time r, and A is the contact area with 
the baseplate. For the contact area we make the 
assumption that the total surface of the aggregate 
interacted after i/ collisions {A = AnE? /v), where 
V is the e-folding width of the exponential function 
in Eq. |4] Indeed, this simple model of maximum 
compaction of a 1/v fraction of the aggregate vol- 
ume only if a previously passive site on the aggre- 
gate surface is hit, yields the very same behavior 
as Eq. |4] which justifies this assumption (see Fig. 
[7]). Thus, we can combine Eqs. [s] and [T7| to calcu- 
late the contact time 



(1-e) 



60 • 2 • TT • i?2 



(18) 



7 




Fig. 7. — A simple 2D sketch of the compacted 
aggregate after some coUisions. According to Eq. 
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39 % of the aggregate radius are compacted 
around an unaltered core. The compacted rim is 
only compressed at the sites where the aggregate 
collided (grey volumes), and a second collision at 
the same site does not lead to further compaction. 
This simple model is capable to explain the in- 
crease of the volume filling factor according to Eq. 
II 

For the parameters of the i? = 1.5 mm aggre- 
gate in section |3.2| this leads to a contact time of 
r = 8.9 ms, which is a realistic result. Indeed, 
the collision time in the experiments of |Heifiel-| 
mann et al. (20071 can be confined to less than 



20 ms (see their image sequence) and preliminary 
studies dedicated to measure the collision time 
of aggregates with (/> « 0.35 with a sohd plate 
yield approximately 5 ms (HeiBelmann et al. , pers. 
comm.). 

We approximate the situation by an elastic 
sphere with a Poisson number of zero, colliding 
with a wall and take the relation for the contact 



time from Hertz (1881) as 



3.85 



(19) 



where G is the shear modulus of the sphere. 



Dintwa et al. ( 2008 1 compile the assumptions 



made in the Hertz model and value the importance 
of frictional contact, non-fiat contact surface and 
large strains. As for the aimed accuracy in our 
model, the deviations from the Hertz model they 
found (and only for large strains) are rather small 
(within a few percent) so that we use Eq 



19 



to 



calculate the shear modulus of the dust aggregate 



to be G = 944 Pa. Combining Eqs. [18] and [19 

we get a scaling relation for the compression of an 
aggregate as 



24.2 • (1 - e) • w • 00 • i?^ 



Vn 



cx K' ■ V 



■0 . „4/5 



(21) 



The upper velocity limit for this extrapolation is 
0.56 m s~^ as for this velocity the volume filling 
factor 0niax = + reaches the physical max- 
imum of (/)RCP — 0.64, which is the random close 
packing of spheres. This high filling factor is, how- 
ever, unlikely to be reached in coUisional compres- 
sion, because the aggregate will then rather frag- 
ment. Thus, the model predicts a fragmentation 
threshold for ~ 0.5 m s~^ in multiple collisions. 
Velocities below a few mm s~^ lead to an insignif- 
icant compaction of A(/)sc ^ 0.01. 

We follow the same Hertzian ansatz to derive a 
scaling relation for the the e-folding collision num- 
ber — AttB^/A. The radius of the contact area 
A in the Hertz model is 



ao = 0.86 • 



G 



(22) 



Thus, the the e-folding collision number i/gc scales 
like 



= 5.42 • 



GR^ 



(23) 
(24) 
M and 



24 



With the scaling relations in Eqs 

/ \4/5 

we get A(psc = • { Q 2m/s J '^sc = 

U ■ ( Q 2ra/s ) ^^'-^ ^^^^ ^'^ S^^*^ ™~ 

crease of the aggregate's volume filling factor in 
each collision as 



A(f>- 



0.2 m/s 



4/5 



\Q.2 m 



-4/5 



(25) 



with 00 = 0.15, A(j) = 0.215, and ly = 850 for 
V < 0.5 m s^^. For a constant velocity this de- 
scription is equivalent to Eq. [4] but it has the ad- 
vantage that it is independent of the collision his- 
tory of an aggregate (e.g. independent of n) and 
is therefore capable to account for non-constant 
bouncing velocities. 
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Table 2: The calculation of collision timescales for different solar-nebula models. For all models we use 
Pd = 300 kg m~^, g = 1 AU, To = 280 K = Mq, and a IQ-^. 

solar-nebula model So [g cm^"^] S R [mm] n [m^"^] v [cm s^^] Tc [years] Tc • i/ [years] 



|Weidenschilling| \l9'77\ 1700 TsO iJo 2.57 • 10"^ 0^27 1^62 1374 

2.66 6.13 10"^ 20.00 0.03 24 

[Andrews and Williams| |2007| 20 0.80 1.50 2.79 • 10"^ 40.20 0.10 85 

0.11 1.79 • 10^ 20.00 0.06 46 

|Desch| 1 (2007) 50500 2.17 1.50 1.40 0.02 4.00 3398 

77.50 7.29 • 10"^ 20.00 0.03 24 



4.2. Collision Timescale 

To value the importance of the bouncing and 
compacting collisions, we estimate the timescale 
on which subsequent collisions occur in the pro- 
toplanetary nebula. For simplicity we make the 
best-case assumption that the entire mass is in 
the relevant aggregate size. A comparable sharp 
size distribution was found for the first growth of 



fractal aggregates ( Krause and Blum 2004J but is 
unknown for the mm-size range. The mean colli- 
sion timescale is then 

Tc= — , (26) 
nva 

where n is the number density of dust aggregates, 
a = AttR^ is the coUisional cross-section of two 
colliding aggregates, and v is the relative velocity. 
A broad size distribution does not extremely alter 
the effect of coUisional compaction. If we consider 
such a wide size distribution and concentrate on 
the compaction of an aggregate at the high-mass 
end by collisions with smaller aggregates, the col- 
lision timescale decreases due to the increasing 
number density of smaller particles, n cx m~^, 
whereas the collision cross section and the rela- 
tive velocity do not significantly change (see Fig. 
[S]). The decrease in collision time is (partly) com- 
pensated by the smaller contact area in the colli- 



sion (see Eq. 22 1 so that the number of collisions 
required to cover the whole surface of the large 
aggregate scales like v (x m~^/^. Thus, the rele- 
vant timescale for the total compaction scales as 
oc m^/'^. Therefore the data given in Table [2] are 
upper limits. If we assume a gas-to-dust ratio of 
100, we can give a general number density of dust 
aggregates in the midplane of the protoplanetary 
disk as 



n = 1.88-10- 



Pdi?3 



1 



Rpda^ 



This equation follows directly from the expres- 
sions for the dust particle density in the midplane, 
Eqs. A16 and A18 inlBrauer et al.|(|2007|). Here, Sq 



is the surface density of the gas in units of [g cm~^] 
at 1 AU in the disk, a is the distance to the star in 
[AU] , S and e are the power indices for the surface 
density and temperature, respectively. Mi, is the 
mass of the star in [Mq], Tq is the temperature at 
1 AU in units of [280 K], and R and are the 
radius in [mm] and mass density in [kg m"'^] of 
a representative dust aggregate. We assume that 
the dust particles are always in equilibrium be- 
tween vertical dust settling towards the midplane 
of the disk and turbulent diffusion which mixes 
the dust up again into the higher regions of the 



protoplanetary disk (DubruUe et al. 1995 Cuzzi 



and WeidenschiUing 2006). Larger particles settle 



closer to the midplane and, hence, lead to higher 
dust number densities. The last square root term 
in Eq. [27| accounts for this effect. 

We assume to be in a nearly non-turbulent re- 
gion in the midplane of the disk, the so-called dead 
zone. Due to the high dust opacity in the midplane 
the ionization degree in this region is too low for 



the magneto rotational instability to operate (Bal- 



bus and Hawley 1991 Reyes-Ruiz 2001). How- 



ever, there are other sources for turbulence, such 



as Kelvin Helmholtz instability (WeidenschiUing 



1979 Johansen et al. 20061, baroclinic instabili- 
ties ( Klahr||2004[ [Petersen et al.||2007[ ), and possi- 
bly free charges mixed to the interior of the disk 
from the upper layers leading to a slight coupling 



of the midplane gas to the magnetic fields ( Turner 
et al.|2007[ ). Therefore, we assume a low, but non- 
zero turbulent a- value of a = 10~^ (Cuzzi and 



Weidenschilling|2006 1 . This low turbulent a- value 



(27) 



influences the number densities of the dust as well 
as the relative velocities of solid particles in the 
disk ( [Ormel and Cuzzi[[2007[ ) . 

We identify three solar nebula models which 
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significantly spread the space of parameters in sur- 
face density. The first model is the minimum 
mass solar nebula (MMSN) model as calculated 
by |Hayashi| ( |1981| ) and |WeidenschilThIg| ( |1977 l. 
We adopt a second model based on recent in- 
terferometric millimeter observation of disks in 
the Taurus-Auriga and Ophiuchus-Scorpius star 



formation regions (Andrews and Williams 20071. 



These observations suggest much flatter surface 
density distributions than in the MMSN model. 
Finally, we consider a revised MMSN model which 
accounts for planetary migration in the early solar 
system ( |Desch|[2007t . In this new MMSN model, 
most of the mass is in the inner parts of the disk 
which leads to very high surface densities of several 
10^ g cm~^, raising the question of gravitational 
instability. The surface densities at 1 AU and the 
power law indices S of each disk model are given in 
Table [S] From these surface densities we can cal- 



culate the number density of aggregates (Eq. 27 1 



also given in Table [2j For this calculation and also 
further on we use Tq = 280 K, e = 0.5, = M©, 
and pd = 300 kg m^'^. 

Different surface densities lead to different rela- 
tive particle velocities in the protostellar disk. Fig- 
ure [8] shows the relative velocities in the midplane 
of the disk at 1 AU for all three models. For cal- 
culating these velocities, we followed [Brauer et al.| 
([2008') and included Brownian motion, relative ra- 
dial drift velocities, and relative velocities caused 
by turbulent gas motion as calculated by |Ormel| 
and Cuzzi (20071. We remark that these relative 



velocities may differ from earlier works due to the 
fact that we adopt more recent calculations of rel- 



ative particle velocities in turbulence (Ormel and 



Cuzzi 20071, and because the solid particle den- 
sity used here is only 300 kg m"'^ accounting for 
porous particle growth. 

To deduce the mean collision velocities of 
two nearly equal sized aggregates with radius 
R, we calculate relative velocities in the inter- 
val accounting for a not perfectly sharp 
size distribution, and take the mean relative ve- 
locity in this interval. Thus, collision velocity, 
number density (Eq. 27 1 and cross section yield 
the collision timescale Tc for different models. We 
also scale the size of the aggregate (cf. section 
4.1) so that we get a mean relative velocity of 
0.2 m s^^. For these sizes and velocities, we get 
collision timcscalcs of less than a year. After the 



time Tc • V, the aggregates are significantly com- 
pacted and all these times are short in terms of 
planet formation. 

4.3. Consequences for Further Protoplan- 
etary Growth 

We address three important consequences of the 
results of this work: (i) The aggregates are com- 
pacted and therefore change their surface-to-mass 
ratio. This has consequences for their coupling to 
the gas and their relative velocities, (ii) The com- 
pacted aggregate possesses macroscopic parame- 
ters like tensile strength, compressive strength and 
shear strength different to the strengths of the 
non-compacted aggregate, (in) The finding of un- 
expected fragmentation requires a review of the 
fragmentation threshold. 

(i) The friction time, the time in which a pro- 
toplanetary dust aggregate is coupled to the sur- 
rounding rarified gas, is in the free molecular flow 
regime deflned as ( Epstein||l924 1 



m 1 



(28) 



where pg and v are the mass density and the mean 
thermal velocity of the gas. m and a a. are proper- 
ties of the dust aggregate, namely, its mass and its 
geometrical cross section (Ja = ttR^ . In section 3.2 



we found that the aggregate volume decreases by a 
factor of two within 1000 collisions without chang- 
ing its mass, which increases its friction time by a 
factor of 1.6. The friction time determines the ag- 
gregate's velocity relative to the nebular gas and, 
thus, relative to other aggregates ( Weidenschilling] 
and Cuzzi]|1993 1 . As the size of the aggregate de- 
creases from the compaction, its relative velocity 
would be that of an uncompacted aggregate with 
twice its diameter. 

(ii) The macroscopic material parameters are 
clearly connected to the coordination number 
(number of contacts per dust grain) and thus 
to the volume flUing factor. The compressive 



strength curve 4>{p) (Blum and Schrapler 2004 



Giittler et al.||2009) paper I) gives the relation for 



the compressive strength as a function of the flUing 
factor. Blum and Schrapler| ( 2004 1 measured the 
tensile strength for differently-compressed dust ag- 
gregates and found an increasing tensile strength, 
closely linear to the coordination number. The 



10 



shear strength (so far not measured for dust ag- 
gregates) is also beheved to be depending on the 



filling factor (Sirono 


2004 Schafer et al. 


2007 


Giittler et al.|2009 


). ( 


jiittler et al. 


(2009 


) perform 



Smooth Particle Hydrodynamics simulations us- 
ing macroscopic material parameters to develop 
a collision model for protoplanetary dust aggre- 



gates. Sirono ( 2004 1 found the occurrence of frag- 
mentation to be depending on the ratio between 
tensile strength and compressive strength. As the 
compressive strength is much more sensitive to 



compaction than the tensile strength (Blum and 
Schrapler||2004 1, the compaction will clearly have 



an influence on the fragmentation threshold which 
is qualitatively shifted to smaller velocities. 

(Hi) The occurrence of fragmentation is rather 
surprising. Earlier experiments (Blum and Miinch 
[1993 , paper II) show fragmentation for veloci- 
ties > 1 m s^^, which is well above the maxi- 
mum velocity of the experiments presented here 
(0.3 m s""'^). One possible explanation is a de- 
creased fragmentation threshold due to the change 
of macroscopic parameters (see (ii)). However, 
p31um and Miinch] ( 1993 ) performed experiments 
with intermediate porosities (0 = 0.26), still with 
a different material (ZrSi04) , and found the same 
threshold. Explanations based on cracking and 
cumulative damage of the aggregate in multiple 
collisions are thinkable to reduce the aggregate 
strength but this - although of major importance 
- remains open for further investigation. A sec- 
ond possibility is a low but non-zero fragmenta- 
tion probability, which would clearly be depend- 
ing on velocity and material parameters, and has 
a finite value Pfrag = 1.8 • 10~^ per collision for 
V « 0.2 m s~^. Although this probability disre- 
gards the history of the aggregate, it is so far the 
only possible treatment of the breakup in multiple 
collisions. 
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turbulence value is a = 10~ 
(dead zone), and, accounting for porosity, the bulk 
density of the aggregates is 300 kg m~^. The ap- 
plied values for the surface density can be found 
in Table E 
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